The plasmon Talbot effect 
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The plasmon analog of the self-imaging Talbot effect is described and theoretically analyzed. 
Rich plasmon carpets containing hot spots are shown to be produced by a row of periodically- 
spaced surface features. A row of holes drilled in a metal film and illuminated from the back side 
is discussed as a realizable implementation of this concept. Self-images of the row are produced, 
separated from the original one by distances up to several hundreds of wavelengths in the examples 

f — , under consideration. The size of the image focal spots is close to half a wavelength and the spot 

positions can be controlled by changing the incidence direction of external illumination, suggesting 
the possibility of using this effect (and its extension to non-periodic surface features) for far-field 

£N| patterning and for long-distance plasmon-based interconnects in plasmonic circuits, energy transfer, 

I and related phenomena. 

Cn INTRODUCTION 

<N 

i_i An important aspect in the development of new branches of optics is the study of analogues of classical optical 

phenomena. In the field of surface waves on metals (surface plasmons) this approach has successfully met with 
engineered nanoscale features providing analogues of lenses and mirrors for future plasmon-based devices [TJ [51 [3] . 
Such manipulation of surface plasmons is of much interest, both from a fundamental viewpoint P~j and from a view to 
applications [21 E]. In practice, it is more difficult to manipulate plasmon fields than their free-space counterparts, as 
they are highly sensitive to metal surface imperfections on the scale of the skin depth (~ 15 nm); nevertheless, they 
have certain advantages, like their ability to concentrate the electromagnetic field near the surface, thus providing a 
I route towards compact light waveguides [1], or their capacity to unveil Raman emission from single molecules through 
enhancement of the local field intensity by several orders of magnitude with respect to the incident light [5]. Here, 
we describe and theoretically analyze the plasmon analogue to another well-known phenomenon of classical optics, 
namely the self-imaging effect discovered by Talbot in 1836 [6] while studying transmission gratings and arrays of holes 

l — 1 perforated in metal films, and later rediscovered and explained by Lord Rayleigh [TJ [H]. The effect is best observed 
through the formation of repeated monochromatic images of a grating at various characteristic distances of the image 
plane with respect to the grating surface. 

Q\ More precisely, a transversally periodic field, paraxially propagating, revives (self-images) to its initial configuration 

after the Talbot distance t = 2a 2 /A, where a is the transverse period and A is the wavelength. In a simple analytical 
description, we represent the grating by a periodic function given in Fourier series form, 

\0 f(x, 0) = f m exp(i2irmx/a), 

o 

where x is the direction of periodicity. The monochromatic wave function emanating from the grating towards the y 
direction reduces then to 

K*^ f(x,y)='Y] fm exp(i2?rmx/a) exp(i27rC m ?//A), (1) 
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where Cm = yl — (m\/a) 2 . The coefficients of x and y in these exponential functions define a vector of magnitude 
27r/ A, the light momentum. In the paraxial approximation (A <C a), the binomial expansion 
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can be truncated at the term proportional to m 2 , equivalent to Fresncl diffraction. This yields 

f(x,y) ps exp(i27rj//A) f m exp(i27rma; / a) exp(-i27rm 2 2//r), (3) 
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from where we immediately deduce 

f(x,r) » exp(i2vry/A) /(x,0), (4) 
f(x,r/2) « exp(i27ry/A) /(a- a/2,0). (5) 

The length r = 2a 2 /A is indeed the Talbot distance at which the initial field self-images (except for an overall phase 
that is washed away when observing intensities), while another image is formed at r/2, laterally shifted by half a 
period and leading to an alternate definition of the Talbot distance [9]. When y is a fraction of r, the field undergoes 
fractional revivals, which in the ideal case are fractal at irrational values of y/r [HI OH HT]. This exotic behavior is 
a consequence of Gauss sums arising from paraxial propagation, which relies on the smallness of the non-paraxiality 
parameter A/a. In practice, this approximation stands only for a finite number of m's in 0, but it can be sufficient 
to render well-defined focal spots, as we shall see below for self-imaging of small features. 

The Talbot effect has been studied in a variety of theoretical and experimental situations [9] [10] HU [12] [13] . 
This phenomenon has an analogue in Schrodingcr evolution of quantum mechanical wavepackets, the quantum and 
fractional revivals of which have been thoroughly discussed [21 [T3]. Although revivals are an exact consequence 
of quantum mechanics, they only arise in optics under the paraxial approximation, and deviation from paraxiality 
destroys the sensitive structure of the Talbot revivals [S]- However, non-paraxial propagation, which only involves 
a finite number of propagating waves, exhibits some good but approximate self-imaging near the paraxial Talbot 
distance [Hfl H7] . 

Self-imaging is not exclusive of periodic objects. The Montgomery effect [TSJIIH] describes for instance perfect image 
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reconstruction of aperiodic objects made of incommensurate harmonic components exp i27r(m/y |m|)x/a 

to replacement of |m| for m 2 in Eq. ([3]), and obviously maintaining the property Q, but not (J5|. Recent work on a 
metal film perforated by quasiperiodic hole arrays has also revealed concentration of transmitted light intensity in hot 
spots at large distances from the film [5D] , suggesting possible extensions of the plasmon Talbot effect to aperiodic 
distributions of surface features. 



SELF-FOCUSING OF PLASMON CARPETS ON METALS: THE PLASMON TALBOT EFFECT 

The analog of the Talbot effect using plasmons is illustrated in Fig. [I] A light plane wave is incident from the back 
of a metal film, planar except for a periodic one-dimensional array of nanoholes or other subwavelength structures, 
with period a. Light is partly transmitted into plasmons on the exit side of the film, thus deploying a complex carpet 
pattern. The field from each of the nanoholes is modeled as a dipolc, oscillating with a frequency corresponding to 
the incident wavelength Ao- This oscillation sets up surface plasmons, propagating into the plasmonic far-field with 
wavelength Asp — Ao/3?{\/e/(l + e)}, which depends on the particular frequency-dependent dielectric function e of 
the metal. We shall concentrate our description on the situation most likely to find practical application, with small 
attenuation and |e| 3> 1, implying that Asp ~ Ao- We shall also concentrate on values of the periodicity a of similar 
lengthscale to the plasmon wavelength Agp. In our graphical illustrations, we model a silver surface with incident 
wavelength Ao = 1.55 fim, for which e = —130.83 + i 3.32 [2T], giving Asp = 1.544 /xm. 

Our detailed analysis starts with the field due to an oscillating single dipole in the y direction at position Ro 
infinitesimally close to the metal surface, incorporating direct propagation and reflection. The electric field, made 
dimensionless through multiplication by Ag, reads [22] 

E s i„gic(r) = J d 2 Q exp [ik • (r - R )] F(Q), (6) 

where 

iA 2 

F(Q) = — ^- [e p k z k y (l-r p )+e s kk x (l + r s )}, (7) 

k = 2tt I Ao is the free-space light momentum, Q = (k x ,k y ) is the projection of the wavevector k into the plane of 
the metal, k z = ^Jk 2 — Q 2 is the component normal to that plane, {k,e p ,e s } is the natural orthonormal basis for 
k, defined as e s = z x k/|z x k| and e p = e s x k, and r p — (ek z — k' z )/(ek z + k' z ) and r s = (k z — k' z )/(k z + k' z ) are 
the appropriate Fresnel reflection coefficients for TM (p) and TE (s) polarization, with k' z = \Jk 2 e — Q 2 [23] . The 
dominant component to E single is E z , and this is strongest on the metal plane for r R in the direction of the dipole 




FIG. 1: Illustration of the plasmon Talbot effect above a metal surface. Light is transmitted through a one-dimensional array 
of nanoholes, setting up a Talbot carpet of interfering plasmon waves. At approximately the Talbot distance r from the array, 
the propagating plasmons revive, giving an array of plasmon focal spots. Plasmon revival at half that distance is also observed, 
with the foci displaced by half the period along the array direction. The dependence of the field on height z above the metal is 
also shown, with the intensity of the z component of the plasmons at fixed height superimposed. The carpet plotted is as for 
Fig.gb). 

(the y direction). Therefore, to maximize the observable effect, we choose to make the periodic dipole array in the x 
direction, with the plasmons propagating in y. 

The ideal plasmon Talbot field comes from an infinite sum of single dipole fields of the form of Eq. ^ , with positions 
at R„ = (na, 0, 0). Using the Poisson sum formula [23], J2 n exp(ik x na) — (2w/a) J2 m &{k x — 2irm/a), the infinite sum 
can be rewritten as a Rayleigh expansion, 

Etotal (r) = — ^2 exp(i2irmx/a) / dk y exp(ik y y + ik z z) F(Q m ) 

m 

= ^ exp(i27rmx / a) F m (y,z), (8) 

n I 

where in the first line Q m = (27rm/a, k y ), and F(Q) is defined in Eq. d7b. In the second line, 2ir/a times the integral 
has been written as the y- and z-dependent Fourier coefficient F m (y,z). Numerical evaluation of this field, for the 
values of the parameters above and various choices of a are shown in Fig. [2ja-c). 

For a — Asp , the Talbot effect is not yet developed, although an interesting periodic pattern appears that could be 
employed to imprint hight-quality 2D arrays. When we move to larger spacing [a — 5Asp in Fig.[2]Jb)], clear evidence of 
self-imaging is observed, which is particularly intense at half the Talbot distance. With even larger spacing [a = 20Agp 
in Fig. |2jc)] a fine Talbot carpet is deployed, showing structures reminiscent of cusp caustics at t and r/2 )25j . The 
focal-spot intensities decrease with distance from the hole array due to plasmon attenuation (as 1.26 mm for silver at 
Ao = 1.55 /mi), to which image contrast is however insensitive at these low-absorption levels. 

The plasmon intensity in the vicinity of slightly less than half the paraxial Talbot distance is shown in Fig. [3] for 
the same conditions as in Fig.^c). The plot on the left shows the field intensity of a focal spot, with cross sectional 
intensities represented on the right. The lateral width of the spot is w 0.5Asp, whereas its extension along y is 
considerably larger. This type of behavior is also observed for other values of the period and for spots at integer 
Talbot distances. The width along x varies from case to case, but it is always close to half a wavelength. 
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FIG. 2: Plasmon Talbot carpets, numerically computed (a-c) from Eqs. |7|-(|8| and analytically approximated (d) from Eqs. 
|7| and (JoJ for different choices of the lattice spacing a: (a) a — Asp; (b) a — 5 Asp; (c,d) a = 20 Asp. The amplitude of the E z 
component of the plasmon field is plotted at a height z = 0.5 fira over a silver surface for a free-space wavelength Aq = 1.55 pm, 
with Asp = 1.544 )im the surface plasmon wavelength. Different scales along horizontal and vertical directions are used in each 
plot: horizontal double arrows show the period a, while vertical arrows signal the paraxial Talbot distance r = 2a 2 /Asp (long 
arrows) and half that distance (short arrows). The hole array is represented by circles in the lower part of each plot. The 
incident light wavevector is along z and its polarization along y (see axes in the center of the figure). 




x=af2-/^. n x=a/2 x=a/2+A co Distance / A. 



FIG. 3: Shape of a plasmon focal spot near half the Talbot distance in Fig. j2|c). The contour plot (left) shows a square of 
side 2Asp centered at (x,y) = (a/2,6), with a — 20Asp and 6 = r/2 — 5Asp = 395Asp. Plasmon intensities at cross sections of 
the spot are given on the right along directions parallel (solid curve) and perpendicular (broken curve) with respect to the hole 
array. 
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FIG. 4: Lattice-period dependence of the intensity near half the Talbot distance at x = a/2 -in the paraxial Talbot effect [Eq. 

the focal spot occurs at exactly y — r/2. The plasmon intensity is represented along y (vertical axis) as a function of lattice 
period a at a height z — 0.5 /im over a silver surface for a free-space wavelength Ao = 1.55 /im. The intensity is normalized to 
the maximum within the plotted range of y for each period. 



ANALYTICAL APPROACH 



It is possible to approximate the field of Eq. (|8| analytically from the observation that the main contribution to 
the integral over k yi particularly in the plasmon far-field, comes from the pole of the r p reflection coefficient, in the Q 
upper-half complex plane. After all, the plasmon dispersion relation derives from that pole (i.e., ek z + k' z = 0), so that 
the plasmon itself posses p symmetry. This contribution may be approximated by the Cauchy integral theorem using 
the r p plasmon pole of wavenumber Qg P = ky/e/(e + 1), corresponding to real plasmon wavelength X SP = 27r/3fJ{Q SP } 

[22) . The remaining z component of the wavevector is fc Zj sp = \/k 2 — Qsp 2 = k/\/l + e. The approximation of the 
integral by the pole residue is appropriate with a cutoff on the Fourier sum in \m\ < N, where N w a/Xsp- For m in 
this range, the approximation gives 

2A (2^) 3 e 2 ( ikz \ ( Ap ^ + T -Cm A fQ . 
F m {y,z) w — exp exp lQspCmy -m , 1 , (9) 

where Cm = \/l — (27rm/aQsp) 2 - This analytical expression yields the same structure as Eq. ([lj), and therefore the 
general explanation of the Talbot effect offered in Sec. applies here as well (assuming that the imaginary part of Qsp 
is small enough to be neglected), apart from the extreme non-paraxiality of the regime under consideration. It should 
be noted that the m dependence of F m is only in the vector and in the exponent of y, and therefore, the Talbot carpet 
is independent of z in this plasmon-pole approximation, except for a global exponential decay away from the surface. 

The evaluation of Eq. ^ corresponding to the conditions of Fig. [2jc) is plotted in Fig. [2jd). Clearly, the ap- 
proximation yields excellent results, particularly in the plasmonic far-field. However, the finite cutoff in the Fourier 
sum implies that there is a finite resolution to all of the interference features in the plasmon field, and hence a finite 
number of fractional revivals (and obviously no fractal revivals), within a Talbot length. 

The choices of the periodicity a < 20Asp in Fig. [2] are in the non-paraxial regime. In Ref. a post-paraxial 
approximation to the classical Talbot effect was studied, in which Eq. ^ was truncated at the term proportional to 
m 4 . The inclusion of this and later terms implies that the field is no longer perfectly periodic, and that the distance 
in y at which the (imperfect) self-imaging occurs is less than t (as in Fig. [2J. However, as our simulations and 
analytic approximation demonstrate, good, if not perfect, Talbot focusing of plasmons should nevertheless be possible 
in practice (similar effects have been noticed in free-space propagation [THE]). The dependence on the period of the 
focal spot near (x,y) = (a/2, r/2), calculated from Eq. is illustrated in Fig. |ij which shows a complex evolution 
of the spot positions, generally below y = r/2. An interesting consequence of these results is that the position of the 
focal spot can be controlled through small changes in wavelength. 

For very large values of |e|, electric dipoles parallel to the metal surface are quenched by their image charges. Then, 
the transmission through the holes depicted in Fig. [TJ relies on parallel magnetic dipoles (provided such dipoles can be 
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induced, for instance under the condition that the metal skin depth is small compared to the hole size [IS]). Magnetic 
dipoles couple best to plasmons propagating in the y direction when they are oriented along x. The above analysis 
remains valid in that case, and in particular Eq. ^ is only corrected by a factor \J e + 1 multiplying the right-hand 
side. Normal electric dipoles (|| z) are also relevant under these conditions, induced by p-polarized light under oblique 
incidence. Again, Eq. (g can be still applied, amended by a factor \fe/Q m . 

DISCUSSION 

Some degree of control over the position of the hot spots is possible when the incident light direction has non-zero 
projection along the hole array direction x: the self-image is displaced along y from the Talbot distance and it is also 
laterally shifted along x, as shown both theoretically and experimentally in Ref. [12]. Under these conditions, the 
projection of the incident light momentum along the hole array, k x , enters Eq. ([9| through an uninteresting overall 

phase factor, but also through the coefficient of the exponential in y, Qsp( m , which becomes y Qsp 2 — (2irm/a + k x ) 2 . 
In the paraxial approximation, one recovers self-imaging at the corrected Talbot distance r x (^sp/Qsp) 3 ; where 

&2/,SP — \l Qsp 2 — (k x ) 2 - Simultaneously, the revival is shifted along x a distance y k x /k y ^p that increases with 
separation from the array. Thus, the position of the focal spots can be controlled through obliquity of the external 
illumination in a setup as in Fig. [I] One should therefore be able to raster the plasmon focus with nanometer accuracy 
for potential applications in nanolithography and biosensing. 

Controllable plasmon focal spots can be particularly advantageous when combined with recently developed adaptive 
ultrafast nano-optics [27], in which femtosecond laser pulses are shaped to achieve a desired objective, such as a time- 
controlled excursion of focal spots in the setup of Fig. [I] 

Superoscillating fields with sub-wavelength localization [28] should also be observed with surface plasmon waves 
using appropriately designed diffraction gratings, as has been recently observed in free-space fields generated by a 
quasi-crystal array of holes [20] . 

The analysis presented here can be straightforwardly extrapolated to other types of 2D light waves, such as guided 
modes in (lossless) dielectric films, long-range surface exciton polaritons |29j . or surface modes in patterned perfect- 
conductor surfaces [30] , with interference between metal patterns and Talbot carpets possibly giving rise to unexpected 
effects in the finer details of the surface modes. The Talbot effect is an attribute of waves, regardless their nature, so 
it must occur in sound, in clastic waves, and in the more exotic scenario offered by electronic surface states in clean 
surfaces like Au(lll), involving wavelengths in the range of a few nanometers at the Fermi level [31] (e.g., Talbot 
carpets could be produced in the vicinity of straight- line steps periodically decorated with adhered nanoparticles). 

CONCLUSION 

We have described theoretically the surface plasmon analogue to the classical Talbot effect. Our numerical calcu- 
lation of the dominant normal component agrees well with our analytic approximation in the plasmon far-field. With 
weak plasmon attenuation, strong focusing of plasmon waves is attainable, even in the non-paraxial regime that we 
have studied, and some control over the position of this focusing is possible by oblique illumination of the incident 
optical field. 

The plasmonic Talbot effect suggests a straightforward and implementable way of tightly focusing plasmon waves 
on a metal surface. Despite the lack of perfect self-imaging imposed by the diffraction limit, the focusing is strong 
enough to allow applications in sensing and imaging. Other potential applications include optical interconnects based 
upon plasmon focal spots aimed at plasmon waveguides. We have emphasized the simplest case in which the effect 
should be strongest, namely the normal component of the field emanating from a periodic array of holes on the metal 
surface. Extensions of the present work to the general case of arbitrary distributions of holes could become an avenue 
to produce on-demand plasmon fields at far distances from the holes. 
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